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Abstract. Meeting time is defined as the time for which two orbits approach each 
other within distance e in phase space. We show that the distribution of the meeting 
time is exponential in (pi, • • • ,pfc)-Bernoulli systems. In the limit of e ^> 0, the 
distribution converges to exp(— ar), where t is the meeting time normalized by the 
average. The exponent is shown to be a = X]/=iP/(l ~ Pi) for the Bernoulli systems. 
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1. Introduction 

In hyperbolic systems, two nearby trajectories separate exponentially in time, which is 
characterized by Lyapunov exponents, the rate of exponential growth of infinitesimal 
initial deviations [H Ej |3]. On the other hand, the Poincare recurrence theorem tells us 
that a trajectory starting at a compact finite size region in phase space returns to the 
initial region infinitely- many times [H |5] . These basic properties characterize complex 
nature of chaotic systems. 

The recurrence time statistics is a reliable tool to measure correlations of chaotic 
trajectories [HI [71 El |9] : if a system is purely chaotic, the recurrence time distribution is 
rigorously shown to be exponential in the limit of small recurrence regions [TOl [TT] . 
In mixed phase space where chaotic components and regular islands coexist, the 
recurrence time obeys a power-law distribution and its exponent is predicted to be 
universal |6l |T2l |T3l |13] , although controversial issues still remain on the existence of the 
universality and the ambiguity of how one determines proper initial recurrence regions. 

In this paper, by analogy with the recurrence time, we will consider time to return 
to not fixed but moving regions. Consider a region which is a neighborhood of a given 
orbit and an orbit starting at the region. One expects that two nearby orbits which 
are initially located at a small but finite-length distance separate and, after a while, 
approach each other within the initial length, possibly infinitely-many times. Statistical 
properties of time intervals for which two orbits come close to each other must contain 
relevant information to chaotic dynamics. In the present study we refer to it as meeting 
time. Although the meeting time is based on the simple idea, to the author's knowledge, 
it is not closely studied as far. In particular, we investigate the meeting time in fully 
chaotic systems. 

In section |2] we define the meeting time distributions for maps with compact phase 
space. Section [3] describes ■ ■ ■ , l/fc)-Bernoulli systems, for which one can obtain 
the meeting time distribution rigorously. In section H] we discuss the meeting time in 
general Bernoulli systems, and concluding remarks are given in section |5l 

2. A definition of the meeting time 

Let X be a phase space and F be a map on X and, for x, x' G X, d{x, x') be a distance. 
Then, for Xi,X2 € X and e > 0, a sequence of times at which the two trajectories 
approach each other within a distance e is given as 




(1) 
(2) 



n« = inf G N I d ^"(xs)) < e, n^~^^ < n} . 



The meeting time is defined as the time interval between these "meetings" : 



(3) 



Then, the meeting time distribution is defined as 



m,(T) = lim ^# {0 < i < iV I T^^ = T} . 



(4) 
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The meeting time depends, in general, on the choice of Xi and X2- If Xi and X2 are both 
periodic orbits, the distribution is trivial. (If simply xi and X2 are both fixed points, 
me(l) = 1, m^(T) = (T > 2) for e > Cc where ec = d{xi,X2) and for e < ec the 
distribution cannot be defined. More generally, if Xi and X2 have different periods, Li 
and L2, Y1t=i'^ ^e{T) = 1, m^{T) = (T > Tl^^l^) for certain e. Here Tl^^l^ is the 
least common multiple of Li and L2. ) 

In ergodic systems, for almost every xi and X2 the distribution is independent of 
initial choice. This is justified as follows: as a function explicitly depending on xi 
and X2, let us here rewrite the distribution as m^{T, xi,X2) and consider a sequence 
of the functions, namely m^lT, F^{xi), F^{x2)) {j = 0, 1, ■ ■ ■). By the definition of the 
distribution, it holds that m^{T, xi, X2) = ■ ■ ■ = m^{T, F^{xi), F^{x2)) = ■ ■ ■■ \yoi every 
T, m^{T, Xi,X2) is an invariant function with respect to Xi and X2 under the map F, 
and, according to the ergodic theorem, is constant for almost every xi and X2- (The 
exceptions are the cases for which xi and/or X2 are periodic.) We hereafter consider 
m^{T) for such xi and X2. 

We will discuss the meeting time by taking the limit e -^^ 0. Since smaller e the 
meeting time becomes longer, the distribution should be taken into account with a 
proper normalization. A plausible and canonical normalization would be an average, 

1 

T, = lim -VtW. (5) 

1=0 

The distribution normalized by the average meeting time is defined as 

m'Xr) = 1# {0 < 2 < AT I tf = t} , (6) 

where t^^ = T^"^ /T,. 

Remark that the meeting time could be regarded as recurrence time, more precisely 
the recurrence for a cross product system of {X,F): let y := X x X be a phase 
space and G be a map on Y, that is G{xi,X2) := (F(xi), F(x2)). The meeting 
time for a system {X, F) is equivalent to the time to return to the diagonal region. 
Re = {{xi,X2) I d{xi,X2) < e}, with respect to the system (Y,G). Thus, infinitely-many 
times meetings are proved to happen due to the Poincare recurrence theorem. 

3. Bernoulli systems 

3.1. (|, -Bernoulli systems 

Let us consider binary symbolic dynamics: let X be a set of semi-infinite binary 
sequences, {0, 1}^, and F be a shift map on X. We here assume (|, ^-Bernoulli 

^This is justified as follows. Here we only show that me{T, xi,X2) = m^iT, F{xi), F(x2)) and 
m^(T,Fi{xi),F^{x2)) = ■m,{T,F^+\xi),F^+Hx2)) follows in the same way. Instead of Eqs. (l)-(3), 
let us here express explicit dependence of n'-*-' and T^'^ on xi and X2 as n'i\xi,X2) and Te^\xi, X2)- 
By the definition, nl^\xi, X2) = ni^\F{xi), F{x2)) — 1 holds and, therefore, the meeting time is given 
as tP{xi,X2) = T^'\f{xi),F{x2)). This leads that m,{T, xi,X2) = m,{T, F{xi), F{x2)). 
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measure. For simplicity we put e = := 2 where r is a positive integer. In 



Appendix A[ we show that the meeting time distribution is given as, 

1/2 (T = l) 

mUT) = { (2 < T < r) , (7) 

4')2-^ (r + l<r) 

(r) 

where J), is an integer generated by the following recurrence relation, 

4^ = j24l, ir + l<T), (8) 

and 4% = 1, J^^ = 2^-'-"^ for r + 2 < T < 2r + 1. One can easily check that 
St=i "^eXT) = 1, and show that the average meeting time is given as 

oo 

T,^ = J2Tm,XT) = r. (9) 

T=l 

The sequence given by the above recurrence relation ([8]) is called the rth generalized 
Fibonacci number. The general solution of the rth generalized Fibonacci number is 
obtained by the roots of its characteristic polynomial 

p^(A) = A" - X'-^ 1. (10) 

Let A(r) denote the largest root of the characteristic polynomial, i.e. Pr{\(r)) = 0. It is 
shown that A(r) is real and its absolute value is greater than 1 and the absolute value of 
the every other roots is smaller than 1 [TS|. Since the polynomial has only one root whose 
absolute value is greater than 1, from general treatment of recurrence relations, one can 
see that the rth generalized Fibonacci number for large T approximately behaves as 

4^ - Af). (11) 
Then, we obtain an asymptotic form of the distribution as 

mUT)^i^^y. (12) 

One can easily show Pr(2) = 1, -Pr(A) > 1(A > 2) and thus 1 < A(r) < 2 [15j. Since 
Pr(2 — 2^^) = 0(2~^), the root A(r.) for large r is approximately given as 

A(,) = 2-2-" + C(2-2"). (13) 

The normalized meeting time distribution is as follows: 

In the limit of large r, namely small e^, the asymptotic form of the limiting distribution 
becomes 

lim m' (r) ~ exp ( ) . (15) 

Figure [T] shows the meeting time distribution obtained by numerical simulations of 
(|, |)-Bernoulli systems. As e decreases, the normalized distribution converges to the 
exponential one of equation (113]) . 
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Figure 1. (a) The meeting time distributions for ^) -Bernoulli system, (b) The 
same plot of (a) with the time normalized as r = T/T^ where the average is 
numerically evaluated. The black line represents the expected one by equation ()15p . 



3.2. {1/k, - ■ ■ , 1/k)- Bernoulli systems 

The derivation made in the previous section can be easily extended to the case of k 
symbols. Let X be {1, ■ ■ ■ , k}^ and assume {l/k, ■ ■ ■ , l//c)-Bernoulh measure and put 
er = k~^. A combinatorial study allows us to derive the meeting time rigorously as well 



as the binary symbol case (see the last paragraph in Appendix A ). The meeting time 
distribution for {1/k, - ■ ■ , l/fc)-Bernoulli systems is given as 

r 1/k (T = i) 

m,,(T) = <^ (2 < T < r) , (16) 

[ jf^^'h-"" (r + l <T) 

where J^^'"-^'' is an integer that satisfies 

r 

jWW = {k-l)J2 4-f (r + 1 < T) . (17) 
i=i 

The average meeting time is given as T^^ = k^. The largest root of the characteristic 
polynomial of the above recurrence relation, denoting it by X(r){k)y becomes, for large r, 
as 

Aw(fc) = k-{k- l)k-'- + 0{k-^'). (18) 

Following the same argument as in the binary case, the limiting normalized distribution 
asymptotically takes the form as 

lim m' {t) ~ exp ( -^^-j-^r] . (19) 

r->oo ^ \ k J 

We should emphasize that the exponent of the limiting exponential distributions 
depends on the number of symbols, while the recurrence time has a simple exponential 
distribution, as far as the system is hyperbolic [TO] . 
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4. General Bernoulli systems 

Next, let us consider {p, g)-Bernoulli systems where 0<p<l, q = 1 — p. Let X be the 
unit interval, X = [0, 1). The (p, g)-Bernoulli map is defined as 

) {x - p)/q {p < X < 1) 

Denote a semi-infinite binary sequence as s = S1S2 ■ ■ ■, where Sj G {0, 1}. The map is 

equivalent to the shift with {p, g)-Bernoulli measure through the following equation, 

00 

(21) 




where /(sj) = <^ * and ips{n) = HiLi fi^i) complement of Sj. In 

<? [Si = i) 

spite of construction of the symbolic expression, the argument in the previous sections 



(mainly developed in [Appendix A] ) cannot be directly applied to this case since the 
probability to appear a given symbol sequence depends on its symbols and p (and q) as 
well. In order to describe the meeting time for {p, g)-Bernoulli systems, we start with 
the case in which xi is periodic and X2 is a non-periodic generic orbit. 

4.1. The meeting time for periodic orbits 

Suppose that xi is a periodic orbit of period L and denote its symbol by t. Obviously, t 
is an L-periodic sequence, i.e. ti = ti+i, and thus it is represented by the first L symbols, 
ti ■ ■ ■ tL. Here, we put e = [/(ti)/(t2) ■ ■ ■ /(^l)]*^- For T > rL + 1, we assume that the 
meeting time distribution for the periodic orbit, denoting it by m^^t (T), satisfies the 
following relation, 

me,t (T + rL) _ ^.,t (T + j) 

ipt{T + rL) ^ MT + J) 

where Q = f{t;)/f{t,). 

The above relation is a generalization of Appendix A to {p, g)-Bernoulli systems. 
(Indeed it is equivalent to equation (|8]) when p = q = 1/2.) It should be noted however 
that the derivation in general cases is not a straightforward extension of (1/2,1/2)- 
Bernoulli systems, but is given via symbolic dynamics with [p, g)-Bernoulli measure 
|16j . We here assume the relation ( 122|) with some arguments which will be more closely 
discussed in section |5l This is justified since xi is periodic and we took appropriate e 
correspondingly. 

Let us consider a polynomial 

r-l 

A''-c^AJ=0, (23) 

j=0 



where c is a constant. In Appendix B , we show that the largest root of the characteristic 
polynomial behaves approximately as Ac = (1 + c) — c(l + c)"'" + O {{1 + c)"^*") for 
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Figure 2. (a) The normalized meeting time distributions with respect to the periodic 
orbit 01011 for the (p, q)-Bernoulli system with p = 0.3. The black line represents the 
expected one by equations and (P7| . (b) The exponent evaluated by fittings as a 
function of p. The lines represent the exponent predicted by (|27l) . Error bars are too 
small to be shown on the plot. 



large r. If the right-hand side prefactor of the recurrence relation f l22|) is constant 
Ci = c, the number given by the recurrence relation at most increases as A^. Since 
we have L-periodic Ci, the asymptotic solution of equation (!22|) for large T is given as 
(T) l^tiT) ~ (Ac^ ■ ■ ■ Xcl)^^^ and then we obtain 



T 



m,,t(r)~(AiA2---AL)^, (24) 

where Aj = 1 — f{ti)f(tiY. Increasing r (namely decreasing e), the factor Aj should 
be normalized by f{tiY, otherwise the distribution diverges. We here assume that by 
replacing Aj by A^ and T by r, the right-hand side of equation ( 1241) gives the 
asymptotic form of the normalized distribution m^^(r). The normalized distribution is 
rewritten as 

1^ rAm)-\..A/fe)-^l^ (25) 



m' 



In the limit of large r (small e), the distribution converges to an exponential function: 
lim t(^) ~ exp f-a^'^V) , (26) 



r— >oo 



where the exponent takes the average over the periodic orbit, namely 

L 



«?'^^ = ^E/(^)- (27) 

i=l 

In figure [2] we show the meeting time distributions obtained by numerical simulations of 
the (p, g)-Bernoulli map, showing good agreement with our theoretical prediction fl26|) . 

In numerical simulations the average meeting time that normalizes the distribution 
is directly evaluated using equation ([5]) and is inversely proportional to e, namely 
Te ~ l/2e. The results imply that the normalization adopted to derive equation (125|) is 
consistent with the numerical results. In our derivations used to provide an asymptotic 
form of the distribution, one can calculate the average meeting time, but one sees that 
this average strongly depends on the parameters (p, q) and the symbols of periodic orbits 
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Figure 3. (a) The normalized meeting time distribution with e — 0.005. (b) The 
decaying exponent of meeting time distributions plotted as a function of p. The curve 
represents the exponent predicted by (1281) . Error bars are too small to be visible on 
the plot. 



(as contrary to simple (l/fc, ■ ■ ■ , 1/ A;)-Bernoulli systems) [IB]. The coincidence between 
the theoretical prediction and numerical observations implies that the normalization 
details are insensitive as far as one takes the limit e — i- 0. 

4-2. The long period limit 

Here we consider the meeting time distribution for which xi is non-periodic. We assume 
that the meeting time for non-periodic orbits is well approximated by the meeting time 
for long periodic orbits. Since the argument in the previous section concluded that the 
average with respect to the period gives the exponent, the distribution for non-periodic 
Xi is exponential for which the exponent is given by the ergodic average. The exponent 
for the long period limit is shown to be 

a^^'''?) = lim a^l'"^ = 2pq. (28) 

Figure |3] shows the distributions obtained numerically for different p. Each exponent 
excellently agrees with the prediction of equation fl28|) when e is sufficiently small. 

4-3. (pi, ■ ■ ■ ,pk) -Bernoulli systems 

The extension of the previous discussions to (pi, • • • ,pfe)-Bernoulli systems is 
straightforward. For Si G §fc := {1, ■ ■ ■ , k}, define a function 

{Pi {Si = 1) 
; , (29) 

Pk {si = k) 

where < < 1 and J2^=iPi — 1- Then, the (pi, ■ ■ ■ ,pfc)-Bernoulli map on the unit 
interval is defined as 

{{x~So)/pi {So<x<Si) 
■. , (30) 

{x-Sk-i)/pk {Sk-i<x<Sk) 
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Figure 4. The exponent of the normahzed meetmg time for (pi,p2,P3)-Bernoulli 
maps. The curves represent the exponent predicted by (|33p . Error bars are too small 
to be shown on the plot. 



where Ss = X]?=iPi- Following the argument in section [3?T| the equation for a periodic 
orbit to satisfy the meeting time distribution is 

m,,t (T + tL) _ (fc) m,^t (T + j) 

^\ '{T + rL) ifl '{T + j) 

where t is the symbol of the periodic orbit and L is its period and y:>i''\n) = nr=i f''''\ti)- 
The constant C^''^ is given as 

For the {pi, ■ ■ ■ , pfc)-Bernoulli map, the long period limit leads to the exponent of the 
normalized distribution as 

k 

aW = 5^pKl-pO- (33) 

1=1 

Note that the obtained exponent is a natural extension of that for {p, g)-Bernoulli 
systems, and even for (l//c, ■ ■ ■ , l//c)-Bernoulli systems. Numerical simulations show 
excellent agreement with our prediction (l33!l . which is shown in figure HI 



5. Conclusions and discussion 



In this paper, we have introduced the meeting time, the time interval for which two orbits 
approach each other within a given distance e. We have shown that the distribution 
of the meeting time is exponential for (l/fc, ■ ■ ■ , 1/ A;)-Bernoulli systems. In the limit of 
small e, the meeting time normalized by its average obeys an exponential distribution 
whose exponent is {k — l)/k. For (pi, ■ ■ ■ ,pfc)-Bernoulli systems, our analysis, based 
on periodic orbit approximations, predicts the exponent as a^''^ = Yl^=iPii^ ~ Pi)- 
This exponent varies in the range of < o;'-'^^ < (k — l)/k and is maximized when 
pi = ■ ■ ■ = Pk = l/k. 
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Figure 5. (a) The meeting time distribution for the piecewise hnear map studied in 
[18] with K = 3. (semi-logarithmic plot). The distribution for every e fits to a simple 
exponential except small deviations due to statistical errors, (b) The meeting time 
distribution for the standard map with K = 1. (double logarithmic plot). We put two 
initial points on the largest chaotic component in the phase space. No quantitative 
changes of the distribution were found even though choosing any other initial points 
in the same component. The distribution obeys a power law partially in particular 
time regimes, while its time scales (and the power law exponents as well) vary as e 
decreases. 



The meeting time distribution is similar to the recurrence time distribution, the 
latter being studied for several chaotic systems not necessarily hyperbolic ones [6l [8]. 
We here point out two aspects on the relation between them. 

First, suppose that Xi is a fixed point, then the meeting time is equivalent to time 
to return to the region of the e-neighborhood of xi, and there is a rigorous proof that 
the recurrence time distribution is a simple exponential [10]. The meeting time contains 
information how xi travels around in phase space. In this sense, the meeting time 
captures more correlations among trajectories than the recurrence time. 

Second, the meeting time distribution for the area-preserving map with mixed 
phase space behaves rather unexpectedly: as shown in figure EJ^a), the meeting time 
distribution in case of the system with sharply-divided phase space is exponential for 
finite- value e. On the other hand, as shown in [171 [IE], the recurrence time distribution 
for the system exhibits power law decay and the corresponding exponent can be derived 
theoretically. Furthermore, the meeting time distribution for a generic system whose 
phase space forms hierarchical mixtures of stable and chaotic regions cannot fit overall 
either to a simple function such as exponential or a power law (see figure EJ^b)). As 
mentioned in introduction, the power law exponent of the recurrence time distribution 
for generic mixed systems is still an unsettled issue, but there exists, at least, a consensus 
that the distribution can be fitted by power law decaying functions in a wide range 
[T2| [T3] . In this way, the meeting time distribution for mixed systems makes a sharp 
contrast, which implies that these two measures capture different aspects of complex 
behavior in hierarchical phase space. 

The origin of the intricate distribution of figure [5](b) is undoubtedly associated with 
complex hierarchical structures. In order to study complex phase space, one may define 
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the meeting time distribution in different ways: for a given point z G X and e > 0, let 
a; be a point for whicli d{z, x) < e. The meeting time for x is defined as 

T,^,{x) ■= m{{nen\d {F''{z), < e} . (34) 

The distribution for z is defined as 

M,,,(T) ■.= fi{{xeX\ diz, x) < e, T,,,(x) = T}) , (35) 

where (■) is an appropriate measure. The defined distribution Mz^eiT) is a point- 
wise measure when one takes the limit e — )■ 0, and thus it essentially depends on z. 
The ergodic average of Mz^^{T) with respect to z is supposed to be equal to m^(T). 
Investigating Mz^e(T) by changing z might be a key to study the case in which phase 
space is inhomogeneous. 

As mentioned in section I4.H the derivation of fl22|) requires the meeting time for 
symbolic dynamics. The meeting time distribution for symbol sequences is defined as 
( l35l) . with the number n of how many symbols coincide, instead of e [H], denoting 
it by Msn{T) where s is a given symbol sequence. In general Ms^n{T) differs from 
Mz,e{T) since the distance for symbolic dynamics is not necessarily equal to that for 
corresponding map systems. If z (and the corresponding symbol sequence t) is periodic, 
the two distribution are identical. In fact, by taking e = [f(ti)f{t2)---f(tL)Y, for a 
point y and the corresponding symbol sequence u, it holds that the first rL symbols of 
u, at least, coincide with those of t iff |?/ — z| < e is satisfied. The same statement holds 
for F\z) (and its symbol expression) with the same e and thus Mz^eiT) = Mt^rLiT). 
What remains to be shown is the following relation 

1 

me,t(r) = -^Mp.(,),,(T). (36) 

1=0 

The distribution in the left-hand side is defined by (jl]) as the long time average of how 
often X2 "meets" one of the periodic points. The right-hand side is an average over 
the distribution with respect to each periodic point. This relation is supposed to be 
justified if X2 uniformly visits the e-neighborhood of one of the periodic points with 
the probability 1/L. These arguments validate that the meeting time distribution for 
periodic points satisfies ( l22l) . 

Recall that the meeting time is the recurrence time for the cross product system, as 
mentioned in the last paragraph of section O While the meeting time is a special case of 
the recurrence time, our numerical observations reveal that the recurrence for the cross 
product system, a couple of two identical dynamical systems, significantly differs from 
the one for the original system. The recurrence time of a typical dynamical system, 
regardless of its dimensionality, has a simple exponential distribution in chaotic systems 
[To] , or a power law one in mixed systems |6]. This in turn suggests that there exist 
dynamical systems whose recurrence time distributions, the recurrence region being 
taken in a special way, do not obey those in generic systems. 

The meeting time that we defined for maps can be similarly considered in 
continuous-time dynamical systems as well. Needless to say, the normalization and/or 
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the small e limit should be properly taken into account in slightly different ways. 
Although we here do not discuss in detail, we expect that the meeting time distribution 
is exponential in fully-chaotic phase space. If this is indeed the case, the exponent may 
have physical significance, or even the relationship to other physical quantities that 
characterize chaotic properties. 

Appendix A. Meeting time distributions for (|, |)-Bernoulli systems 

Let X be {0, 1}^ and denote a semi-infinite symbol sequence, s e X, as 



where Sj e {0, 1}, and F be the shift map, i.e. F{s) — S2S3S4 ■ ■ ■. For s, s' e X, the 
distance is defined as 



Putting e = er := where r is a positive integer, one can see that s meets s', namely 
d{s, s') < Cr, iff the first r symbols of s and s' coincide, i.e. Si = s'^,i = 1, ■ ■ ■ , r. 

Taking account of the fact that every symbol, Si or s-, is either '0' or '1' with the 
probability 1/2, the meeting time is given as follows: let s be fixed without loss of 
generality and the first r symbols of s' are those of s and for i > r -|- 1 is either Si or 
where Si is the complement of Sj. Let us consider a symbol sequence, t^+i ■ ■ ■ tr+n, where 
ti is either Si or 'si. For a given t^+i ■ ■ ■ tr+n, let [tr+i ■ ■ ■ tr+n] denote the subset of s' such 
that s^_,_j = tr+i, i = 1, ■ ■ ■ , n and call it as a sequence of length n. The probability to 
appear a symbol sequence in [tr+i ■ ■ ■ tr+n] is 2~", denoting it by ([tr+i ■ ■ ■ tr+n]) = 2~". 

The meeting time distribution m^^{T) is given by the probability of symbol 
sequences of s' such that, after T-times shift, the first r symbols coincide with those of 
s, namely the r symbols after the Tth symbol coincide = st+i, i — 1, • • • , r. For 
T — 1 the distribution is given by [s^+i] i.e. me^{l) — ^{[sr+i]) — 1/2. For T > 2 a 
symbol is a subset of [s^+i]. For 2 < T < r, me^{T) — since Sr+i is in the r symbols to 
coincide. For T > r + 1, the distribution me^{T) is given by finite sequences of length 
T with following two properties: 

(i) Any of the last r symbols are not the complement (the last r symbols are expressed 

as [• • -St+I • --ST+r])- 

(ii) There are no r-consecutive symbols of the non-complement expression before the 
last r symbol sequence, (the r non-complement symbols appear only in the last of 
the whole sequence, otherwise its meeting time is smaller than T.) 

Let here denote the number of such sequences of length T. One finds that J^^J^ = 1 
since [sr+iSr+2 ■ ■ ■ S2r+i] is the unique symbol sequence for T = r-|- 1, and J^^^2 = since 
[sr+iSr+2Sr+3 ' ' ' S2r+2] for T = r -|- 2. For T > r + 3, consider the following sequence. 



S1S2S3 • • • , 



(A.l) 



00 




(A.2) 



i=l 



Sr+1 tr+2 ■ ■ ■ tr-l StSt+1 ' ' ' ST+r- 



(A.3) 
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For T < 2r + 1, the bracket part tr+2 ■ ■ ■ ^r-i allows complete binary combinations since 
its length is less than r, and therefore j!p^ = 2^^'"^^. 

(r) 

The number sequence Jj. , which is the combinatorial number of tr+2- ■ - tr-i 
containing no r-consecutive sequences of s7, satisfies the recurrence relation given as 

J? = E 41, (r + 1 < T) . (A.4) 
This is simply justified as follows: let us consider a set of the binary sequences {0, 1} 

(r) 

of length i which do not contain r-consecutive 'O's. Denote this set by and 
nf''^ = ^Bf\ (n-^^ is obviously equivalent to 4^ by replacing with '0' and Si '1'.) 
In order to construct B\ , let us start with the simplest case r = 2: the set of length 

(1'] f2) il') 

1 is i?{ = {0, 1}. B\ is generated from -Bj„\ by the following two rules: (a) append 
'1' to every element of Bf^-^ and (b) append '0' to an element of Bf\ if the last symbol 
of the element is not '0'. Then, Bl is recursively generated, e.g. i?2 = {01, 11, 10}, 
-83^-' = {Oil, 111, 101, 010, 110}, ■ ■ ■. The appending rules (a) and (b) lead the recurrence 

relation, nf^ = nf^-^ + ^j-^2- obviously n^-j*^ represents the number of elements generated 

(2) (2) 
by (a). The number of elements by (b) is ni_2 since the elements of Bl_\ whose last 

symbol is not '0' are generated from elements of by (a). For r > 3, B^^^^ is generated 

from -Bj-I.^ as well by (a) and the modified rule of (b): (b') append '0' to an element of 

Bl_\ if the last r — 1 symbol sequence is not (r — l)-consecutive 'O's. (e.g. for r = 3 

= {0,1}, = {01,11,00,10}, fif^ = {011,111,001,101,010,110,100}, ■■■.) 

The number of elements for which the modified rule (b') is applied is n'f_}2 + ■ ■ ■ + nf^r 

since an element of B^^^ generated from bI%, ■ ■ ■ , B^l^^. by (a) contains at least one 

'I's in the last r — 1 symbol sequence, that is indeed the requirement of the modified 

rule (b'). Thus, one obtains the recurrence relation nf^ = Y7j=i ^t-j ^^"^ equation (IA.4p . 

Before closing the Appendix, we make a remark on the case of k symbols. The 
settings for fc-symbol cases are simply given, as described in section 13.21 Since the 
complement of Sj has k — 1 candidates of the symbols, the recurrence relation which the 
numbers of the sequences that hold (i) and (ii) in fc-symbol cases satisfy has the factor 
— 1 in the right-hand side of equation (]A.4p . (In other words, the rule (a) in /c-symbol 
cases allows to append any of {1, 2, ■ ■ ■ , A; — 1} so that the factor k — 1 appears.) 

Appendix B. The real largest root of equation ( I23|) 

Let us rewrite the left-hand side of equation ( 12^ as 

P,,,(A) = A'- - c^l^, (B.l) 

where c > is a real constant and r > 1 is an integer. One can easily show Pr.,c(l + c) = 1 
and for > 0, 

[i + c + dY -I 



Pr,c{l + c + d) = {l + c + dy -c- 



c + d 
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= [(1 + c + dy - 1]-^ + 1 > 1. (B.2) 

One finds that the real largest root of the polynomial is smaller than 1 + c. 

The derivative of Pr-,c(A) at A = 1 + c is P;^^(l + c) = ^^^f^ > and is of order 
of (1 + cy. Therefore, for large r, the largest root is supposed to be close to 1 + c and 
can be obtained by perturbative calculations. Let us rewrite the polynomial, by the 
expansion around 1 + c, as 

P,,e(l + c + x) = P,,e(l + c) + P;,,(l + c)x + ^P;;,(l + c)a;' + • • • 

+ ipW(l + cK, (B.3) 
and expand x by the power series of (1 + c)~^ as 

x = ai(l + c)-^ + 02(1 + c)-^' + ■■■. (B.4) 
Then, the expansion of the polynomial becomes as 

P,,,(l + c + x) = 1 + il±^-ll[ai(l + c)-' + 02(1 + c)-2^ + •••], 

+ l^ci^ + + c)-^ + a2(l + c)-'' + . . .]2 + . . . , 

= l + ^ + C»((l + c)-'-). (B.5) 

The root is obtained by calculating Oj such that every order of (1 + c)"^ becomes zero. 
One obtains Oi = — c for the lowest order and the root as (1 + c) — c(l + c)~^ H . 
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